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Generic Kerr orbits exhibit intricate three-dimensional motion. We offer a classification scheme 
for these intricate orbits in terms of periodic orbits. The crucial insight is that for a given effective 
angular momentum L and angle of inclination t, there exists a discrete set of orbits that are geomet- 
rically n-leaf clovers in a precessing orbital plane. When viewed in the full three dimensions, these 
orbits are periodic in r — ^. Each n-leaf clover is associated with a rational number, 1 -I- qre — uje /uJr, 
that measures the degree of perihelion precession in the precessing orbital plane. The rational num- 
ber qrB varies monotonically with the orbital energy and with the orbital eccentricity. Since any 
bound orbit can be approximated as near one of these periodic n-leaf clovers, this special set offers 
a skeleton that illuminates the structure of all bound Kerr orbits, in or out of the equatorial plane. 
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I. INTRODUCTION 

Black hole orbits are defined by precession. The per- 
fectly closed ellipse of Kepler's Laws gives way to the rel- 
ativistic precession of Mercury's perihelion in the weak 
field around a star. In the strong-field, perihelion preces- 
sion in the equatorial plane of a black hole can result in 
zoom-whirl orbits for which the precession is so great at 
closest approach that the particle executes multiple cir- 
cles before falling out to apastron again. An orbit out of 
the equatorial plane, the plane perpendicular to the spin 
axis of the black hole, is shaped by yet another kind of 
precession - precession of the orbital plane. These most 
general black hole orbits live in three dimensions, are not 
confined to a stationary plane, and are dynamically intri- 
cate. A complete classification of these 3-D orbits is the 
purview of this article. 

Carter famously showed that there were four constants 
of motion[TJ[2] for the orbits of spinning black holes, one 
for each canonical momentum, so that the orbits are inte- 
grable. Still, black hole orbits have long evaded a simple 
geometric classification. While any geodesic orbit could 
be computed easily, a concise general account of how 
changes to the constants of motion would alter its shape 
was unavailable. Recently a topological taxonomy based 
on periodic orbits provided a complete classification of 
all equatorial orbits [3]. 

In brief, Ref. [3] shows that just as Mercury is a pre- 
cession of the ellipse, any relativistic orbit can be under- 
stood as a precession of a periodic orbit. Although there 
is no ellipse in relativity, no 1-leaf clover, there are 2- 
leaf, 3-leaf,... n-leaf clovers as well as n-leaf clovers with 
nearly circular whirls. The equatorial periodic orbits are 
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defined by a rational number 

<Z.^ = ^ - 1 (1) 

where is an average angular frequency in the equa- 
torial plane and ujr is the radial frequency. Aperiodic 
orbits correspond to irrational ratios of frequency while 
periodic orbits correspond to rational qr^p. The number 
qrip explicitly measures the degree of perihelion preces- 
sion beyond the ellipse as well as the topology of the 
orbit. The q^p = 1/3 orbit is a 3-leaf clover while the 
Qrip = 1 + 1/3 orbit is a 3-leaf clover with 1 whirl per 
radial cycle. And, importantly, the qrip — 1/3 -I- e orbit 
looks like a 3-leaf clover precessing at a rate of 27re of 
azimuth per radial cycle. (For a complete description see 
Ref. [Sill].) The classification is especially effective since 
qrip varies monotonically with the energy of an orbit for 
a given L. As the value of qrip increases, the topology of 
the orbit varies in a systematic way as the energy and 
orbital eccentricity also increase (for a given L). The re- 
sulting taxonomy nicely exposes the complete equatorial 
dynamics. 

The goal here is to generalize the equatorial taxonomy 
to fully generic 3-D Kerr motion. We could identify fully 
periodic orbits and argue that all generic orbits are ap- 
proximated at arbitrary precision by that set of measure 
zero [3]. However, it is sufficient to consider the less re- 
strictive, larger set of orbits that are periodic only in 
as these will be shown to be perfectly periodic when pro- 
jected into an instantaneous orbital plane, as illustrated 
in Fig. [1] 

A series of orbits is shown in 3-D on the leftmost col- 
umn of Fig. [TJ in the r — cos 9 plane in the middle col- 
umn, and projected in an effective orbital plane in the 
final column. These orbits are closed in r — 6 and also 
in the orbital plane, but are not fully closed in 3-D. The 
following sections will be devoted to realizing this argu- 
ment. Similar reasoning led to a taxonomy of generic 3-D 
orbits in a Post-Newtonian expansion of two black holes 
in Ref. [SIE]. 
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FIG. 1: r-6 periodic orbits with L = 3, E = 0.932516, cost = 0.4 and o = 0.99, but with different r — phasing. Column 1 
shows the full 3D orbit. Column 2 is the projection of the orbit into the r-cosO plane. Column 3 is the projection into the 
orbital plane. All rows have ro = Va ~ 8.82713 and ipo = 0. The initial 6 values are as follows: Row 1 9o = Omin = 0.414139; 
Row 2 00 = 0.8; Row 3 = f ; Row 4 So = 2; Row 5 9o = 6lmax = 2.72745. 
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II. THE BASICS 

We begin with the Kerr metric in Boycr-Lindquist co- 
ordinates and geometrized units (G = c = 1) and the con- 
ventional choice oi M — jj. ~ 1: 



(2) 

-dtdif +^dr^ + me"^ 



2r\ ,.9 4arsin^0 ,. , ^ ■ 2 , i/i2 



A 



9 „ / 9 9 2a^r sin^ \ , 9 
sure H + ' ' "'■"'^ 



dip'^ 



where 



E = + cos^ 6 
A = - 2r + 



(3) 



Carter reduced the equations to first integrals of motion 
[1112], exploiting the four constants of motion E,Lz,Q 
and /i: 
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We will often refer to equations ([4) as the Carter equa- 
tions. In those equations, an overdot denotes differenti- 
ation with respect to Mino time [Tj, A, which is related 
to the particle's proper time, r, by d\ 
quantities 



and the 
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(5) 
(6) 



are the polar and radial quasi-potentials, respectively |8]. 

The quasi-potentials reveal some well-known geomet- 
ric information about bound non-plunging orbits (orbits 
that neither escape to infinity nor cross the horizon of the 
central black hole). First, they reveal the radial turning 
points, which occur at roots of R{r). For a given E,Lz 
and Q, the quartic polynomial has four roots. The outer- 
most two are periastron and apastron, between which the 
radial position of a bound orbit oscillates. A similar anal- 
ysis of the roots of 8(6') reveals that every bound orbit os- 
cillates between a fixed ^min and 6'inax symmetrically dis- 
tributed about the equatorial plane^, i.e. ^min = ti" — ^max 
[8HTO]. The upshot is that every 31? orbit will gener- 
ally lie in a toroidal wedge around the equatorial plane 



^ For equatorial orbits, 



bounded Tp and in radial coordinate and between 6'niax 
and TT — flmax in polar angle [11] . 

Every bound Kerr orbit also has an associated triplet 
of fundamental frequencies (wj., "^e, w,p) , which can be 
defined for any choice of time coordinate [T^]. The sim- 
plicity afforded by the choice of Mino time and exploited 
heavily in [T] [TH] is that, since the radial and polar mo- 
tions decouple in Mino time, each of r(A) and 9{\) is inde- 
pendently periodic. As a result, the Mino-time frequen- 
cies can be defined and computed directly from equations 
dial and (|4b|. 



We will only be concerned with the radial and polar 
frequencies here. To obtain them, we first define the 
radial and polar Mino periods via 
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The corresponding Mino-time frequencies are then 
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(7b) 
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Note that we use Mino time purely for ease and conve- 
nience and that the frequency ratios which figure promi- 
nently in our analysis are independent of the choice of 
time variable. 

We want to consider 3D orbits that are closed in r — 
9. That closure will result when the polar and radial 
frequencies are rationally related, or in language more 
directly useful for our orbital plane description of the 
motion, when the quantity 



uje_ 

UJr 



1 



(9) 



= 7r/2. 



is rational. To be useful, a classification based on orbits 
with rational qre has to have two properties: the rational 
Qre must tell us about the topology of the orbit, and 
it must relate that topology to more physical conserved 
quantities. In the subsequent sections, we show that this 
is indeed the case. 



A. The Energy Spectrum 

In the spirit of the equatorial classification of we 
begin by describing how q^e varies with energy. For ease, 
and in anticipation of the fact that our analysis will ul- 
timately focus on the discrete set of E values for those 
orbits with rational values of (7^9, we will refer (loosely) to 
the dependence of E on g^e as an energy spectrum. The 
subtlety in establishing a simple relationship between qre 
and E is the choice of which other parameters to keep 
fixed as E is varied. In the appendix we show that the 
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FIG. 2: Left: Five VcB curves with different L values but all with a = 0.99 and cost = 0.4. In order from the highest to 
the lowest curve, the L values are: L = 3.4 > Libso, L — Libso = 3.32432, Lisao < L = 3 < Libao, L = iisso = 2.85501, and 
L = 2.7 < Lisso. The horizontal line at VcS — \ shows the energy of marginoUy bound orbits. Right: Five VcS curves with 
different L values but all with a = 0.99 and cosi = —0.4. The highest curve has an L = 4.5 > Libso- The second highest curve 
has L — Libso = 4.28330. The middle curve has Lisso < L = 4 < Libso- The second shortest curve has an L = Lisso ~ 3.74594. 
And the smallest curve has L = 3-5 < Ligso- The horizontal line at Ves = | shows the energy of the marginoUy bound orbits- 



key combinations are an effective total angular momen- 
tum L and angle of inclination t for orbits around a black 
hole of a given spin a defined by 

L2 = L2 



Q 



(10) 



cos I : 



L 



first used by [131 US] and used occasionally in other ref- 
erences [ini El HH HZ] - our construction turns out to be 
greatly facilitated by varying L while keeping t fixed, as 
opposed to varying while keeping Q fixed. 

This choice of orbital parameters allows us to write 
equations ([s]) and (|6| as 

Q{0) = L^sm\-cos'9\a\l-E^) ' 



sin^6i 
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(12) 

With this particular combination of constants, we can 
produce an analog of the familiar Schwarzschild effec- 
tive potential for nonequatorial Kerr motion. Consider a 
black hole of given spin a. For a non-spinning black hole 
(a = 0), we can rewrite the radial equation (4a) as 

2 



1 fdrV 



(13) 



This standard effective potential formulation of 
Schwarzschild motion relates the radial velocity with 
respect to particle proper time to an effective energy 
Ecff — E^/2 and an effective potential T4ff that is a 
different function of r for each fixed L - crucially, VcS 
is independent of E. The result is a simple visual way 
to describe the different types of allowed motion as L is 
varied. 

However, for fully 31? orbits around a spinning black 
hole, an analogous potential is not self-evident. The 
counterpart to equation ( |13| (which must involve veloc- 
ities with respect to Mino-time in order to decouple the 
radial motion from the polar motion) is 



1 



- ? = 



(14) 



In this case, the dependence on E in R{r) cannot be sim- 
ply separated and moved to the right-hand side. It would 
seem that the best we can do with eqn. ( 14 ) is end up 



with a Veff that depends on all the constants of motion. 
We therefore lose the ability to visualize easily the varia- 
tion of orbits with energy, as we can in the Schwarzschild 
case, because even at fixed Lz and Q (or L and i), chang- 
ing E also causes the potential to shift. As written, then, 
equation ( 14 ) admits a one-dimensional effective poten- 



tial description, but that description is not useful because 
there is a different potential for every combination of or- 
bital parameters. 

However, if we consider only the behavior at the turn- 
ing points, we can construct a useful pseudo-effective po- 
tential. The idea is to set r = in Eqn. (14), which 
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amounts to setting R{r) = 0, and to solve for E: 



E (a, i, L, r) 



2a cos lLt + {a? + (-2 + r) r) (r^ {L^ + r^) + a2 (2 + r) {L^ - cos^ lL^ + r^)) 

r{r^ + a2 (2 + r)) 



(15) 



We then define a pseudo-effective potential 

KffUo = |^ ■ (16) 

that will allow us to draw various lines of fixed _B on a 
potential that maintains its shape and visually identify 
turning points of the motion, even if the difference be- 
tween E and the value of VcS no longer gives the value 
of r^. 

Fig. [2] illustrates the utility of this approach. For every 
fixed L, as L is lowered we get an analogous pattern of 
orbits to the Schwarzschild case. There is a minimum of 
Vee, which corresponds to a stable constant radius orbit. 
There is a maximum of Vcff, which corresponds to an 
unstable constant radius orbit. Unlike the Schwarzschild 
case, the constant radius orbits are not circles. Instead 
they lie on the surface of a sphere bounded between 0,nax 
and ^inin = TT — 6'max- Wc hcrcaftcr call these spherical 
orbits. 

The parallel story continues. As L is lowered, the un- 
stable spherical orbit becomes bound once a certain crit- 
ical L value is crossed. That value Libso is the angular 
momentum of the innermost bound spherical orbit (ibso) , 
the unstable spherical orbit with critical energy E — 1. 
An innermost stable spherical orbit (isso) appears as a 
saddle point of Vcs once L drops to yet another critical 
value -Lisso- For L < Ljsso, all orbits plunge into the cen- 
tral black hole. Fig. [2] demonstrates the consistency for 
both prograde and retrograde orbits. 

If we had chosen to keep Lz fixed while varying Q 
instead of keeping L fixed while varying t, we would not 
have seen the same simple pattern. Appendix [B] shows 
the breakdown in the Schwarzschild analogy when using 
orbital parameters {Lz,Q). 

The result, for a given L, t, is that q^e increases mono- 
tonically with energy. The lowest energy bound orbit is 
the stable spherical orbit, and, importantly, this orbit 
has the lowest value of qre for that combination of L, t. 
As detailed in Ref. [3], the constant radius orbits do not 
have rational value zero, as can be proven by taking the 
zero eccentricity limit, e — > 0. 

Since qre is monotonic, its upper bound q^'^^ is the 
value of qre for the maximum energy bound non-plunging 
orbit for a given L. Whether q^^^ is finite or infinite 
depends on whether L is greater than or less than Libso- 
If L > Libso, the unstable spherical orbit is unbound and 
has energy E > 1. q'^^^ is therefore the qre value of the 
E = 1 orbit, and despite the fact that the E = 1 orbit just 



reaches r = 00 after infinite time, its qre is nonetheless 
finite. As we reduce L, q'^g^^ increases monotonically, and 
eventually q^g'^^ — )• 00 once L = Libso- For all L < Libso, 
q'^e^^ remains infinite 121 US] • This happens because the 
maximum energy bound non-plunging orbit is now the 
homoclinic orbit (or separatrix orbit), which formally has 
an infinite number of whirls during its lone infinite-period 
radial cycle. A detailed analysis of the homoclinic orbit 
can be found in P| HH] • 

Figure |3] is a plot of the qre versus energy for a given 
a, L and 3 sets of L values. It is representative of the 
general trend we see for any combination^ of a,L,L. As 
the energy increases, so does qre- As L decreases towards 
iisso, the minium value of qre increases. This trend was 
seen equatorially in Ref. ^3J. 

In figure [3] we see that the qre also increases with ec- 
centricity, e. Again this is a general trend so that qre is 
monotonic with eccentricity. The larger qre, again for a 
fixed (a,L, l), the more eccentric the orbit. 

We have shown that qre corresponds to an energy spec- 
trum for 3-D orbits. What we want now is to show this 
also corresponds to a measure of zoom-whirliness and so 
is also a toplogical indicator. As we will see, quite incred- 
ibly, this qre measures the amount by which the angle in 
the orbital plane overshoots 27r, that is, precesses, in one 
radial period. In other words, when qre is rational, it is a 
direct measure of the topology of the orbit in the orbital 
plane and increases monotonically with energy, thereby 
defining a spectrum of zoom-whirl orbits in the orbital 
plane. 



B. Periodic Tables and the Orbital Plane 

We preface this section with the caveat that the orbital 
plane construction below naively employs flat space vec- 
tor algebra and vector calculus constructions (e.g. cross 
products of 3-vectors) without fully taking into account 
the curvature of the background Kerr spacetime. Prima 
facie, it is not obvious that the formalism should accu- 
rately capture geometric or topological features of 3D 
orbits. Nevertheless, we have the amazing result that 
the r — 6 periodic orbits correspond to a spectrum of 



^ The case of t = 0, tt needs to be handled as in Ref. [3] because 
that is motion that takes place entirely in the equatorial plane. 
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FIG. 3: Top: The plot shows the monotonic relationship be- 
tween Qre and energy for all bound orbits with a given a, L 
and cos t. We show three different L values all with a = 0.99 
and cos t = 0.4. The graphs cut off on the left at the en- 
ergy value for the stable spherical orbit with that a, t and 
L. Bottom: The plot shows, for the above parameter values, 
the monotonic relationship between qre and orbital eccentric- 
ity e = y^j^- The lower eccentricity bound is e = 0, also 
corresponding to the stable spherical orbits. 



zoom-whirl orbits in this effective orbital plane, beau- 
tifully mirroring the equatorial result of Ref. [3]. For 
now, we simply state our results, which are compelling, 
and report that a more precise analysis of the connection 
between the orbital plane construction and a relativisti- 
cally precise projection of the motion using local tetrads 



is underway. A very precise implementation for the PN- 
expansion of two black holes can be found in Refs. [3 [6]. 

We consider the projection oi r — 9 periodic orbits in 
an instantaneous orbital plane that we define naively as 
the plane in the tangent space spanned by R and P, 
defined below, with a corresponding angular momentum 
C = R X P. At every instant, the orbital plane is the 
plane perpendicular to the angular momentum vector. 

It is useful to define 

p^ir^ + a^/^ (17) 

and convert from ellipsoidal to Cartesian coordinates 

X =p sin cos (fi 
y =p sin 9 sin </? 
z =r cos 9 . 



Then, 
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where i = x,y,z and k,j = r,9,ip. For convenience we 
take the M ^ limit HE, 
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(22) 



To find the orbital plane, we write 
so that we can define 

X^kxL 
Y =ixX . 



(23) 



(24) 
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The orbital plane is spanned by X, Y. (For a more de- 
tailed exposition on the orbital plane variables, see Ref. 
[SI [3].) This informally defined orbital plane is sufficient, 
as we will see, since it effectively soaks out any ip motion. 

Fig. [4] shows a table of orbits in the effective orbital 
plane. Our periodic table assembles orbits with rational 
Qrff as an energy spectrum, with energy increasing from 
top to bottom and then from left to right. The topol- 
ogy of zoom-whirl orbits in the effective orbital plane is 
encoded in qre through 

V 

q^g + - , (25) 

where w is the number of nearly circular whirls and v 
indicates the order in which the z zooms, or leaves, are 
traced out. So the qro = 1 + 2/3 orbit is a (z = 3)-leaf 
clover, that executes w = 1 whirls during each each radial 
cycle before it moves to the v = 2 leaf in the pattern. 

This result is quite remarkable: q^e is a measure of the 
number of times the orbit returns to ^min per radial cy- 
cle, yet it gives topological information about the degree 
of precession in a very different angular variable, namely 
the angle swept out in the oribtal plane. Had we instead 
projected the orbit onto the r — cos 9 plane, our r — 9 pe- 
riodic orbits would look like Lissajous figures as in Fig. [5] 
The geometric information in Fig. |4] is severely obscured 
when the trajectories are plotted as Lissajous figures. 

Fig. [T] shows trajectories with the same orbital param- 
eters but different r — 9 phasing. All orbits have the 
same E, L, l and therefore the same (r^, Tp, 0max)- How- 
ever, Ta coincides with different initial values of 9 in the 
range tt — 6'inax < 9o < ^max for each picture. Under 
shifts in r — phase, the iD orbits are all rather different 
(illustrated in the first column) as are their correspond- 
ing Lissajous figures (illustrated in the second column). 
Notice, in stark contrast, that varying the initial phas- 
ing of r-vs.-9 merely corresponds to an overall rotation 
of the very same zoom-whirl orbit in the orbital plane 
(illustrated in the final column). 



Moreover, as Fig. [T] illustrates, a change in r — 6* phase 
corresponds to a simple rotation of the orbit in the effec- 
tive orbital plane. An orbit that hits apastron at 6'„iax 
will be rotated by tt/2 in the orbital plane relative to 
an orbit with identical (E, L, u) that hits apastron at 
9 = 7r/2. 

Any aperiodic orbit will be arbitrarily well- 
approximated by a nearby periodic orbit. What's 
more, aperiodic orbits will look like precessions of 
low-leaf clovers. Just as Mercury is a precession of the 
ellipse, an orbit with q^g = 1/2 -I- e is the precession of a 
2-leaf clover that accumulates an extra 27re of azimuth 
during each radial cycle. Our results therefore provide a 
complete taxonomy for generic Kerr orbits. 
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III. SUMMARY 



Our results are neatly summarized in Figures [4] and [T] 
Fig. |4] illustrates that orbits periodic in r — 6* assemble 
into a spectrum of multi-leaf clovers when projected in 
a loosely defined orbital plane. The topology of the or- 
bit is encoded in a rational number q^-g = — — 1, from 
which one can immediately read off the number of leaves 
(or zooms), the ordering of the leaves, and the number of 
whirls. For a given L, t, the rational number q^e mono- 
tonically increases with energy and with eccentricity. So, 
a simple 3- leaf clover [q^g = 1/3) has less energy and is 
less eccentric than a 2-leaf {q^-g = 1/2) of the same L,t. 
Significantly, the rational number q^e is bounded below 
so that there are no qre — >■ orbits in the strong-field 
regime. There are therefore no tightly precessing ellipti- 
cal orbits in the strong-field regime. All eccentric orbits 
will have a countable number of leaves. 
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Appendix A: Spherical Orbits 

In the cases of Schwarzschild and equatorial Kerr mo- 
tion, orbits of constant r — circular orbits — serve to or- 
ganize the ranges of orbital parameters over which bound, 
nonplunging motion exists. Constant r orbits in the gen- 
eral Kerr geometry play a similar organizational role but 
need not lie in a plane. Thus, they are not necessarily 
circular orbits but rather spherical. Spherical orbits were 
first treated in [S] and later analyzed in the context of ra- 
diation reaction in [10l[T6] (in the latter references, these 
constant r orbits are refered to as "circular, nonequa- 
torial orbits" , but we use the original shorter moniker 
"spherical" from Ref. [S]). Like circular orbits, spherical 
orbits have f = f = 0; unlike their circular counterparts, 
spherical orbits do not have ^ = 0. 

An initial analysis of equatorial Kerr motion (we can 
think of Schwarzschild motion as the a — Q subcase) be- 
gins with expressions for and E of circular orbits as 
a function of r and the (fixed) central black hole spin 
a. Our generic Kerr analysis will reproduce one such 
equatorial-like picture for each inclination l and will have 
an effective total angular momentum L take the place of 
the more conventional conserved quantity but oth- 
erwise proceed analogously. We therefore turn now to 
deriving expressions for the effective angular momentum 
L and E of spherical orbits as functions of r, a and l. 
Ref. [TD] has similar expressions for Q and Lz of spher- 
ical orbits in terms of r, a and E^ but as we explain in 
Appendix |B] aggregating orbits with fixed values of the 
constants t and L is most conducive to a clear exposition 
of the dynamics. 

As in the equatorial Kerr case, our starting point is 
the radial quasi-potential R{r). We begin by expressing 
R{r) and its derivatives in terms of E, l and L. From 
eqn. pi. 



The condition f = implies R{r) =0 from equation 



( 4a I . Solving for f from equation ( 4a I we find that 



R (r) = {E^ - 1) r* + 2r^ + (a^ [e^ - l] - L^) 

+ 2r {a^E^ - 2aEL cos l + L^) + a^L^ (cos^ l - l) 

(Al) 

R' (r) ^4{E^ - 1) + 6r^ + 2 (a^ {E^ - l} - L^) r 
+ 2 (a'^E'^ - 2aEL cos t + L'^) 

(A2) 

R" (r) = 12 {E^ - 1) + 12r -|- 2 (a^ {E^ - l} - L^) ■ 

(A3) 



1 R 

1 fR' 
271 

2^ 



(A4) 



where R' {r) ^ We can see immediately from equa- 
tion ( [Ai] ) that f = implies R' (r) = 0. Similarly, 'r = 
implies that R" (r) = 0. 



To find expressions for all Eg and Ls for a fixed a and 
L, we set R (r) = R' (r) = and solve for Eg (r, a, t) and 
Ls {r,a,i). Solving the two coupled quadratic equations 
yields four solutions for each of Eg and Lg. We determine 
the physically admissible solutions by imposing that Lg 
always be positive, i.e. an effective angular momentum 
magnitude. Additionally, because each fixed l should 
replicate the orbital structure of the Schwarzschild ge- 
ometry, both the Lg and Eg solutions should asymptote 
at low r- values to the innermost time- like spherical orbit. 
There should also be a minimum Lg and Eg value corre- 
sponding to the innermost bound spherical orbit (ibso). 
And the r at which the minima occur on the Lg and Eg 
graphs should be the same. Finally, at large r, our Lg 
plot should reproduce the Newtonian limit, ^/L oc r and 
Eg should asymptote to 1. 



Combining the above conditions, we find 
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Es [r, a, i) 



(-3 + r) (-2 + rf + a* sin* t (1 + r) 



- 2ar cos tA (-a^ sin^ l + r^) y r (-a"* sin^ l + 2a? sin^ lA + r^) 

- a^r^ sin^ t[a^ {4-4 (-1 + r) r + cos^ i (1 + r) (-5 + 4r)} 
+ 2 (-1 + r) r {2 - 3 (-2 + r) r + cos^ t (-4 + r (-1 + 2r))} ] 



+ aV^ [4 (-2 + r) {1 + (-3 + r) r} + cos^ t{8 + r (-23 + (17 - 4r) r) }] 



(A5) 



-a^ sin^ t - 2a2r2 sin^ i - r'^) x 



I 



-^Jr (-a^ sin^ t + 20^ sin^ tA + r'^) + ar cos t (a^ + r (-4 + 3r)) 

Lg (r, a, l) — ^ Eg (r, a, t) 

-a4 sin^ t - (-2 + r)^ r2 + (4 - 2r + cos2 t (-3 + 2r)) 



|- (-3 + rf r^ - sin^ t (1 + rf + 20^^ (- (-3 + r) (1 + r) + cos^ t (-3 + r^)) | 



(A6) 



We recover the functions Ec and ic given in [20] for equa- 
torial Kerr circular orbits by setting t = for prograde 
and t = TT for retrograde in equations (A5| and (A6). 
From there, we recover the well-known Schwarzschild 
functions Ec an d Lc (see, for i nstance, Ref. [19]) by set- 
ting a = in (A5) and (A6) (note that, by spherical 
symmetry, those values must be and are independent of 
t). 

Figure |6] shows both Lg and Eg as functions of r with 
parameters cost = 0.4 and a = 0.99. The following qual- 
itative features are representative of all t and a values 
and mimic the features of Schwarzschild. Both Eg and Lg 
have minima that occur at the same r. The minimum Lg, 
^issoj corresponds to the least Lg for which there exists a 
spherical orbit. The V^g plot corresponding to L = Lisso 
has a saddle point where the stable and unstable spheri- 
cal orbits merge. For all L > Ljsso there are two spheri- 
cal orbits, whose r-values exactly correspond to the local 
minimum and maximum of the effective potential plots 
of that L, t and a. The maximum is the unstable spher- 
ical orbit and the minimum is the stable spherical orbit. 
There is a critical value Lg ~ iibso at which the unstable 
spherical orbit has Eg — 1, and for all Lg > Libso, the 
unstable spherical orbit is unbound with Eg > 1. For a 
fixed t and a, all the qualitative properties of the generic 
Kerr orbits replicate the Schwarzschild system. 

The innermost bound spherical orbit, ibso, is defined 
as the spherical orbit with critical energy i?ibso = 1. To 
find the Libso and nbso, we set ( Al ) and ( A2 1 to zero with 
E = 1. The innermost stable spherical orbit, isso, is the 
minimum of the Lg plot and is subject to the further 
constraint R" (r) — 0. We therefore find the isso for a 
given t and a by setting all three of equations ( Al ), ( A2 ) 



and ( A3 ) to zero simultaneously and solving for L;, 



Appendix B: Choosing conserved quantities 

The Kerr metric has four conserved quantities. They 
are conventionally chosen to be the black hole mass (/i), 
the orbital energy (i?), the z-component of angular mo- 
mentum (Lz) and the carter constant (Q). Because each 
of those quantities are constants of the motion, any com- 
bination of them is also a constant of the motion. There- 
fore, there are an infinite number of choices of four in- 
dependent quantities we could make for our conserved 
quantities. 

We have chosen to use ^, E, effective angular momen- 
tum (L, where L = ^/Q~+~L^) and inclination angle (t, 



where cos 



). This section provides an explanation 



for our choice. 

Our goal was to realize a generic Kerr orbit struc- 
ture that generalized the Schwarzschild and equatorial 
Kerr orbit structures presented in [3]. To bring that 
goal to fruition, we look for a set of conserved quanti- 
ties such that we could hold one fixed and reproduce all 
the qualitative features of Schwarzschild dynamics (isso, 
ibso, etc.). 

Using the conventional Q, L^ and E, the equatorial 
Kerr system is defined by Q = 0. There are two sets 
of Lz and E solutions for circular orbits, one prograde 
and one retrograde. Figure [7] shows the two solutions for 
a — 0.995 and Q — Q. We can see that the solutions never 
intersect and each solution has all the qualitative features 
present in the standard organization of Schwarzschild or- 
bits. 

However, when Q becomes large enough, regardless of 
the spin, we see a loss of adherence to these features. 
Specifically, there is no longer an isso, and the two sets 
of solutions for Lz and E for a fixed Q mix. While this 
phenomenon is not seen until Q gets large, it is present 



12 



L' ' I - for all spin values. The discontinuity in the and E 

spherical graphs, as well as the loss of the isso is seen for 
8 - - the full range of a values. 



3 - 



2' 
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r 
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0.94 - 



0.92 - 




0.90 I 
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r 

FIG. 6: Top: The figure shows a plot of Ls vs r for spherical 
Kerr orbits with a = 0.99 and cos l = 0.4. Bottom: Shows 
a plot of Es vs r for spherical Kerr orbits with the a = 0.99 
and cos l — 0.4. 



The upshot is that there are values of Q that do not 
allow us to reproduce the familiar qualitative organiza- 
tion of Schwarzschild dynamics if we choose to look at 
orbits of constant Q as an ensemble. In contrast, we 
find that with {E,L, t), for every fixed t, the qualitative 
dynamical picture mimics the familiar Schwarzschild one 
beautifully. In this picture, each t corresponds to a fixed 
orbital inclination so that equatorial orbits correspond to 
one of two t values: l = for prograde, and l = tt for ret- 
rograde. Furthermore, whereas each fixed Q admits two 
associated E and solutions each for spherical orbits, 
each L produces only one curve each for Lg and Eg. 

Figure [8] shows a set of E and Lz plots for spherical 
orbits with Q — 12.5. We can see the loss of the isso and 
the mixing of the two seperate solutions. The curves are 
no longer even single-valued at a given r. Moreover, the 
Es{r) curve can have more than 2 orbits with a given 
E, as opposed to only the stable and unstable constant 
r orbits we are used to in the Schwarzschild effective po- 
tential picture. We have picked four points on the fixed 
Q plots, each with a unique set of orbital parameters, E, 
Lz and Q. For each of those points, we have determined 
the corresponding E, L and l and plotted the Eg[r) and 
Ls{r) curves for each of those t values. Notice that there 
is no such breakdown when we look at curves of fixed 
t rather than fixed Q. Instead, the latter curves faith- 
fully reproduce the expected qualitative features of the 
corresponding Schwarzschild or equatorial Kerr curves. 
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FIG. 7: Left: The figure shows a plot of Lz vs r for equatorial circular Kerr orbits with a = 0.995 and Q = 0. Right: Shows a 
plot of S vs r for circular Kerr equatorial orbits with a = 0.995 and Q = 0. 
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a=0.995, Q=12.5 
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FIG. 8: The above pictures show the organizational difTerences between using the conserved quantities E, and Q and using 
E, L and l. AU plots arc for spherical orbits with a — 0.995. Top: Curves of L^^ vs. r and Eg vs. r for spherical orbits all with 
fixed Q = 12.5. Below: Curves of Ls{r) and Es{r) for spherical orbits with the four different fixed t values associated with the 
four labeled points. Each such point corresponds to the same orbit in both the fixed Q and fixed t graphs. Their parameter 
values are: (1) E = 0.99, = 0.598971, r = 3.01492, L = 3.58591, i = 1.40298; (2) E = 0.98, = -2.28682, r = 5.09346, 
L = 4.21065, L = 2.14493; (3) E = 0.955, = -1.49806, r = 8.92632, L = 3.83981, t = 1.97158; (4) E = 0.97, = 2.47180, 
r = 15.9948, L = 4.31391, i = 0.960654. Orbits 1 & 2 are unstable; 3 & 4 are stable. 



